Geodesic flows on semidirect-product Lie groups: geometry of singular measure-valued solutions

The EPDiff equation (or dispersionless Camassa-Holm equation in 1D) is a well known example of geodesic motion on the Diff group of smooth invertible maps (diffeomorphisms). Its recent two-component extension governs geodesic motion on the semidirect product ${\rm Diff}\circledS{\cal F}$, where $\mathcal{F}$ denotes the space of scalar functions. This paper generalizes the second construction to consider geodesic motion on ${\rm Diff} \circledS\mathfrak{g}$, where $\mathfrak{g}$ denotes the space of scalar functions that take values on a certain Lie algebra (for example, $\mathfrak{g}=\mathcal{F}\otimes\mathfrak{so}(3)$). Measure-valued delta-like solutions are shown to be momentum maps possessing a dual pair structure, thereby extending previous results for the EPDiff equation. The collective Hamiltonians are shown to fit into the Kaluza-Klein theory of particles in a Yang-Mills field and these formulations are shown to apply also at the continuum PDE level. In the continuum description, the Kaluza-Klein approach produces the Kelvin circulation theorem.Comment: 22 pages, 2 figures. Submitted to Proc. R. Soc.

GG-Strands

A $G$-strand is a map $g(t,{s}):\,\mathbb{R}\times\mathbb{R}\to G$ for a Lie group $G$ that follows from Hamilton's principle for a certain class of $G$-invariant Lagrangians. The SO(3)-strand is the $G$-strand version of the rigid body equation and it may be regarded physically as a continuous spin chain. Here, $SO(3)_K$-strand dynamics for ellipsoidal rotations is derived as an Euler-Poincar\'e system for a certain class of variations and recast as a Lie-Poisson system for coadjoint flow with the same Hamiltonian structure as for a perfect complex fluid. For a special Hamiltonian, the $SO(3)_K$-strand is mapped into a completely integrable generalization of the classical chiral model for the SO(3)-strand. Analogous results are obtained for the $Sp(2)$-strand. The $Sp(2)$-strand is the $G$-strand version of the $Sp(2)$ Bloch-Iserles ordinary differential equation, whose solutions exhibit dynamical sorting. Numerical solutions show nonlinear interactions of coherent wave-like solutions in both cases. ${\rm Diff}(\mathbb{R})$-strand equations on the diffeomorphism group $G={\rm Diff}(\mathbb{R})$ are also introduced and shown to admit solutions with singular support (e.g., peakons).Comment: 35 pages, 5 figures, 3rd version. To appear in J Nonlin Sc

Commutator errors in large-eddy simulation

Commutator errors arise in large-eddy simulation of incompressible turbulent flow from the application of non-uniform filters to the continuity -- and Navier-Stokes equations. For non-uniform, high order filters with bounded moments the magnitude of the commutator errors is shown to be of the same order as that of the turbulent stress fluxes. Consequently, one cannot reduce the size of the commutator errors independently of the turbulent stress terms by any judicious construction of such filter operators. Independent control over the commutator errors compared to the turbulent stress fluxes can, instead, be obtained by appropriately restricting the spatial variations of the filter-width and filter-skewness. For situations in which the dynamical consequences of the commutator errors are significant, e.g., near solid boundaries, explicit similarity modelling for the commutator errors is proposed, including the application of Leray regularization. The performance of this commutator error parametrization is illustrated for the one-dimensional Burgers equation. The Leray approach is found to capture the filtered flow with higher accuracy than conventional similarity modelling, which is particularly relevant for large filter-width variations

Formation and Evolution of Singularities in Anisotropic Geometric Continua

Evolutionary PDEs for geometric order parameters that admit propagating singular solutions are introduced and discussed. These singular solutions arise as a result of the competition between nonlinear and nonlocal processes in various familiar vector spaces. Several examples are given. The motivating example is the directed self assembly of a large number of particles for technological purposes such as nano-science processes, in which the particle interactions are anisotropic. This application leads to the derivation and analysis of gradient flow equations on Lie algebras. The Riemann structure of these gradient flow equations is also discussed.Comment: 38 pages, 4 figures. Physica D, submitte

Alpha model for 3D Eulerian mean fluid circulation

We provide a summary of the known analytical properties of the Alpha models, including an outline of their derivation and the associated assumptions, their simplification for the case of constant dispersion length (alpha) and their conservation properties. We also offer interpretations of nonlinear dynamics of the viscous alpha models and indicate the differences one might expect from the dynamics of the Navier-Stokes equations

Lyapunov stability of ideal compressible and incompressible fluid equilibria in three dimensions

Linearized stability of ideal compressible and incompressible fluid equilibria in three dimensions is analyzed using Lyapunov's direct method. An action principle is given for the Eulerian and Lagrangian fluid descriptions and the family of constants of motion due to symmetry under fluid-particle relabelling is derived in the form of Ertel's theorem for each description. In an augmented Euleriah description, the steady equilibrium flows of these two fluids theories are identified as critical points of the conserved Lyapunov functionals defined by the sum, H + C, of the energy H, and the Ertel constants of motion, C. It turns out that unconditional linear Lyapunov stability of these flows in the norm provided by the second variation of H + C is precluded by vortex-particle stretching, even for otherwise shear-stable flows. Conditional Lyapunov stability of these flows is discussed. 24 refs

A minimal no-radiation approximation to Einstein's field equations

An approximation to Einstein's field equations in Arnowitt-Deser-Misner (ADM) canonical formalism is presented which corresponds to the magneto-hydrodynamics (MHD) approximation in electrodynamics. It results in coupled elliptic equations which represent the maximum of elliptic-type structure of Einstein's theory and naturally generalizes previous conformal-flat truncations of the theory. The Hamiltonian, in this approximation, is identical with the non-dissipative part of the Einsteinian one through the third post-Newtonian order. The proposed scheme, where stationary spacetimes are exactly reproduced, should be useful to construct {\em realistic} initial data for general relativistic simulations as well as to model astrophysical scenarios, where gravitational radiation reaction can be neglected.Comment: 9 page

Energy efficiency of integrated crop management systems in the Dark Brown Soil Zone

Non-Peer Reviewe

On asymptotically equivalent shallow water wave equations

The integrable 3rd-order Korteweg-de Vries (KdV) equation emerges uniquely at linear order in the asymptotic expansion for unidirectional shallow water waves. However, at quadratic order, this asymptotic expansion produces an entire {\it family} of shallow water wave equations that are asymptotically equivalent to each other, under a group of nonlinear, nonlocal, normal-form transformations introduced by Kodama in combination with the application of the Helmholtz-operator. These Kodama-Helmholtz transformations are used to present connections between shallow water waves, the integrable 5th-order Korteweg-de Vries equation, and a generalization of the Camassa-Holm (CH) equation that contains an additional integrable case. The dispersion relation of the full water wave problem and any equation in this family agree to 5th order. The travelling wave solutions of the CH equation are shown to agree to 5th order with the exact solution

Analytic solutions and Singularity formation for the Peakon b--Family equations

Using the Abstract Cauchy-Kowalewski Theorem we prove that the $b$-family equation admits, locally in time, a unique analytic solution. Moreover, if the initial data is real analytic and it belongs to $H^s$ with $s > 3/2$, and the momentum density $u_0 - u_{0,{xx}}$ does not change sign, we prove that the solution stays analytic globally in time, for $b\geq 1$. Using pseudospectral numerical methods, we study, also, the singularity formation for the $b$-family equations with the singularity tracking method. This method allows us to follow the process of the singularity formation in the complex plane as the singularity approaches the real axis, estimating the rate of decay of the Fourier spectrum